Dynamics of composite Haldane spin chains in IPA-C11CI3 
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Magnetic excitations in the quasi-one-dimensional antiferromagnet IPA-CuCl3 are studied by cold 
neutron inelastic scattering. Strongly dispersive gap excitations are observed. Contrary to previously 
proposed models, the system is best described as an asymmetric quantum spin ladder. The observed 
spectrum is interpreted in terms of "composite" Haldane spin chains. The key difference from actual 
S — 1 chains is a sharp cutoff of the single-magnon spectrum at a certain critical wave vector. 



Antiferromagnetic 2-leg spin (5 = 1/2) ladders and 
the closely related 5 = 1 Haldane spin chains Q, are 
an example of quantum disorder and mass generation in 
extended spin networks. In these quantum spin liquids, 
the magnetism is suppressed due to collective zero-point 
fluctuations and the unique topology of one dimension 
(ID), rather than to finite system size. Spin ladders 
play a key role in the dynamics of stripe phases in high- 
temperature superconductors 3] and under certain con- 
ditions can themselves support exotic types of supercon- 
ductivity |4| . On a more fundamental level they are ideal 
models for studying the collective spin dynamics in ID, 
quantum critical points and phase transitions in external 
magnetic fields, and the effects of quenched disorder. 

Despite the wealth of relevant theoretical results, ex- 
perimental studies are lagging behind for shortage of suit- 
able model compounds. In the best known ladder ex- 
amples Sri4Cu2404i 5] and SrCu203 6] the large en- 
ergy scales of magnetic interactions limit spectroscopic 
studies, especially at high fields. A typical problem with 
many known Haldane-gap systems [3, |S| is a large single- 
ion magnetic anisotropy that is often associated with 
5=1 spins, and qualitatively affects the dynamics and 
field behavior 9]. Moreover, a Haldane spin chain in- 
trinsically has a smaller Hilbert space than 5 = 1/2 
spin ladders, and therefore lacks certain very interest- 
ing spectral features 10]. In the present work we re- 
port the discovery of an isotropic ladder spin network 
in the 5 = 1/2 compound (CHa^CHNHsCuCls (IPA- 
CUCI3 ), that was previously thought to be a prototyp- 
ical ferromagnetic-antiferromagnetic (F-AF) spin chain 
[lH . We use inelastic neutron scattering to study its 
magnetic excitation spectrum. The experiments reveal a 
spectacular truncation of the magnon branch at certain 
critical wave vectors that we attributed to peculiarities 
of magnon interactions, and bring new insight on diluted 
(CH 3 )2CHNH 3 Cu(Cl a .Bri_ a; )3 compounds [H . 

IPA-CUCI3 crystallizes in a triclinic space group PI 
with a = 7.766 A, b = 9.705 A, c = 6.083 A, a = 97.62°, 
(3 = 101.05°, and 7 = 67.28° 11]. The key features of 
the structure are shown in Fig. ^ The magnetism is due 
to 5 = 1/2-carrying Cu 2+ ions arranged in sheets par- 



FIG. 1: Layers of magnetic Cu 2+ (red) and the bridging CP 
ions (cyan) in IPA-CUCI3 . Previously proposed alternating 
Cu chains (dotted line) run along the c axis. Actual asym- 
metric spin ladders (shaded) are parallel to the a direction. 
The ladder rungs are defined by the vector d. 



allel to the (a, c) crystallographic plane. These sheets 
are well separated by non- magnetic organic layers. A 
spin gap of A w 1.5 meV in IPA-CUCI3 was first dis- 
covered by Manaka et al. [ll| in x(T) measurements. 
This value is consistent with the critical field H c ~ 11 T 
that induces an ordered antiferromagnetic phase at low 
temperatures 13]. The gap was attributed to a singlet 
ground state of bond-alternating Cu 2+ chains running 
along the crystallographic c axis, as shown by dotted 
line in Fig. Indeed, this model consistently explained 
all magnetization curves and ESR experiments, with a 
quantitative agreement obtained assuming alternating F- 
AF bonds ^]|. An F coupling between nearest-neighbor 
Cu 2+ sites Ji < 14] is consistent with the relevant Cu- 
Cl-Cu bond angles in the crystal structure of IPA-CUCI3 . 

At the center of the present work are inelastic neutron 
scattering experiments on deuterated IPA-CuCk sam- 
ples, prepared by crystallization from solution 11]. The 
x(T) and ESR data for the deuterated samples agree well 
with those for the non-deuterated ones. 20 single crystals 
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FIG. 2: Time of flight (a) and 3-axis (b)-(d) inelastic neutron data measured in IPA-CuCl3 at T = 1.5 K. In (a) the suppression 
of inelastic intensity at h < 0.5 in (a) is due to a large / (momentum transfer along e*) and the FM-dimer structure factor for 
the magnon branch (see text). For h > 0.5, / is small at all times : |/| < 0.2. In (b) - (d) the shaded areas are calculated peak 
shapes due to resolution. 



with total mass of 3.5g were co- aligned to a mosaic spread 
of 4° . The data were collected using cold-neutron 3-axis 
spectrometers NG5-SPINS in NIST and IN12 in ILL, and 
at the Disc Chopper time-of-flight Spectrometer (DCS) 
at NIST 15]. The sample was mounted with the (a, c) 
crystallographic plane parallel to the scattering plane of 
the instruments and maintained at or below T = 1.5 K. 
On SPINS we utilized (guide) - 80' - 80' - (open) colli- 
mations and a BeO filter after the sample for a fixed final 
neutron energy Ef = 3.7 meV, selected by a flat pyrolitic 
graphite (PG) analyzer. On IN-12 a Soller collimator of 
60' was used only at pre-sample position. A Be filter was 
positioned after the sample. Neutrons of fixed- incident 
E[ = 4.7 meV or fixed- final Ef = 3.5 meV energies were 
used in conjunction with a horizontally focusing PG ana- 
lyzer. A constant background (typically 1.5 counts/min, 
depending on configuration) was subtracted from all 3- 
axis data sets. On the DCS instrument the data were 
taken with E\ = 6.7 meV neutrons, the incident beam 
forming an angle of 60° with the a* direction. The back- 
ground was directly measured by removing the sample 
from the cryostat. 

A series of constant-^ and constant-g scans revealed 
well-defined long-lived magnetic gap excitations in a large 
part of the Brillouin zone in IPA-CUCI3. Typical data 
are shown in Fig. Our measurements unambiguously 
show that the magnetic strong-coupling direction in IPA- 
CUCI3 is along the crystallographic a axis. The global 
dispersion minimum is at q = ( 2n 2 fl , 0, 0) with n-integer 
(Fig. EK) 5 and also corresponds to a maximum of ob- 
served inelastic intensity. In contrast, the dispersion 
along the c* axis is rather weak, with an energy min- 



imum at Z = in Fig. Ob- The originally proposed 
model of bond-alternating F-AF chains running along 
the crystallographic c axis is thus totally inconsistent 
with our results. Nevertheless, the measured gap energy 
A ~ 1.2 meV is in agreement with that deduced from 
bulk measurements. 

In the wave vector range 0.2 < h < 0.8 and 1.2 < h < 
1.8, where sharp magnon peaks were observed, the 3- 
axis data were analyzed using a model single-mode (SM) 
cross section, with an empirical spin-gap dispersion re- 
lation 0, slightly modified to accomodate a transverse 
disperion: 

S(q,w) oc S SM (q)(^-^ q ), (1) 
uo\ = Uq cos 2 (tt/i) + [A 2 + 4bl sin 2 (7r/)] sin 2 (7r/i) + 
+ c 2 sin 2 (27r/i), (2) 

where q = /ia* + /cb*+/c*. The gap is equal to A at / = 0, 
and is allowed to slightly oscillate along c* with a trans- 
verse bandwidth bo. The boundary energy corresponds 
to coo. The cross section was numerically convoluted 
with the spectrometer resolution function calculated in 
the Popovici approximation [l7^ and fit to the experimen- 
tal data. The magnetic form factor for Cu 2+ ions was 
built into the fits. Separate coefficients 5' SM (q) were used 
for each constant-q scan to determine the 3-dimensional 
structure factor in the most model-independent man- 
ner possible. The parameters of the dispersion relation 
in Eq. Q were treated as global for all the collected 
scans. The fitting procedure yielded c^o = 4.08 (9) meV, 
A = 1.17 (1) meV, b = 0.67 (1) meV and c = 2.15 (9) 
meV. Scans simulated using these parameter values are 
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FIG. 3: Magnon dispersion measured in IPA-CuCl3 along the 
a* (a) and c* (b) directions (symbols). Heavy solid lines are 
fits, as described in the text. In (a) the green shaded areas 
indicate the domain of the 2-magnon continuum. Solid arrows 
show the positions of scans in Fig. El (c) and (d). Dashed 
arrows show the critical wave vectors. 

shown as shaded areas in Figs.O>d. The dispersion rela- 
tions are plotted in solid lines in Fig. 03 Symbols in these 
plots were obtained in fits to individual scans, as opposed 
to global fits to the measured data. Figure^ shows the 
/-dependence of the single-mode component of the dy- 
namic structure factor 5 SM (q) measured at h = 0.5. The 
/i-dependence of this structure factor for Z = is shown 
in Fig. QJ}. 

The measured c-axis modulation of the magnon inten- 
sity is exactly explained by a ferromagnetic dimer with 
a particular spin separation d in the crystal structure 
(Fig.DJ: d= (-0.0854, -0.1316, 0.4432). The data on 
the h = 2n 2 fl reciprocal-space rods are very well repro- 
duced by assuming 5 SM (q) oc cos 2 {nhd x + 7rld z ) / 'a; q , as 
shown by the solid line in Fig.^K- To a good approxima- 
tion one can thus describe IPA-CUCI3 in terms of uniform 
chains made up of these ferromagnetically correlated spin 
pairs. The chains run along the crystallographic a axis 
and are composed of effective 5=1 spins. Since the 
single mode intensity has its maximum at h = 2n ^~ 1 , 
(Fig. Et>), the correlations between the nearest-neighbor 
effective 5 = 1 spins are predominantly antiferromag- 
netic. If one assumes that near the magnetic zone-center 
a dominant fraction of the total spectral weight is con- 
tained in the single mode excitations, one can repro- 
duce the measured /^-dependence of intensity with only 
nearest-neighbor AF interactions in the 5=1 chains. 
For this model S(h) oc sin 2 (71-/1)/^, which is in good 



agreement with experiment (solid line in Fig. 0Jd). We 
conclude that the singlet ground state in IPA-CUCI3 and 
the gap excitations are due to composite 5 = 1 Haldane 
spin chains that run perpendicular to the originally pro- 
posed bond- alternating chains directions. 

As shown in Fig. the crystal structure of IPA- 
CUCI3 presents a number of possible AF superexchgange 
routes between Cu 2+ spins along the a axis. In the most 
general case, the appropriate model for IPA-CUCI3 is a 
spin ladder with F rung interactions Ji, a leg coupling 
constant J2 and a diagonal exchange interaction of mag- 
nitude J3 (presumably also AF). If this interaction was 
infinitely strong, the mapping onto a Haldane spin chain 
would be exact. In fact, the measured dispersion differs 
from that for an ideal 5 = 1 chain. In the latter, the 
ratio of the spin wave velocity v to the gap energy A 
is v/A ~ 6 [13 . In contrast, in IPA-CuCl 3 v/A « 3.1, 
about twice as small. The analogy with Haldane spin 
chains is thus only qualitative, and the |Ji| is probably 
of the same order of magnitude as J 2 and J3. Without 
a direct comparison with first-principle calculations, the 
exact values of the exchange constants can not be de- 
duced from the measured dispersion and intensity mod- 
ulation of the single-magnon excitations. The observed 
small dispersion of magnetic excitations along the c* di- 
rection indicates a weak coupling J4 between individual 
ladders (Fig.EJ. From the position of the dispersion min- 
imum we conclude that J4 < 0. 

A remarkable feature of the measured spectrum is the 
abrupt disappearance of the magnon branch at a cer- 
tain critical wave vector q c . This phenomenon is well 
illustrated by the time of flight data of Fig. EK The 
sharp gap excitation extends to h = | + h c ~ 0.8, where 
h c ~ 0.3, and then vanishes abruptly, rather than per- 
sisting to the zone-boundary at h = 1. Constant-g 3-axis 
data reveal similar behavior at the equivalent wave vec- 
tor /i = | — /i c ~ 1.2. The well-defined resolution limited 
peak seen at h = 1.2 (Fig. Eh) is weakened and broadened 
at h = 1.15 (Fig. EH, where the shaded area is the simu- 
lated profile for a resolution-limited peak), and is totally 
absent at h = 1.1 (not shown). 

Such a dramatic spectrum termination phenomenon 
has been recently discovered and investigated in the 
quasi-2D material PHCC 19]. To our knowledge, IPA- 
CUCI3 is the only other known spin gap material with a 
similar feature. The effect can be explained by a cou- 
pling between single-magnon and multi-magnon states 
|lfj |2£| . The domain of the two-magnon continuum in 
IPA-CUCI3 , though invisible in our experiments due to 
low intensity and a very tight energy resolution, can be 
calculated using the measured one-magnon dispersion. 
As represented by the green shaded areas in Fig. Si, the 
continuum has a gap of 2 A. Note that the single-magnon 
branch crosses into the 2-magnon continuum at exactly 
h = 2n + 1 ± h c . Energy and momentum conservation 
laws forbid any single-magnon excitation outside the do- 
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FIG. 4: Measured /- (a) and h- dependencies (b) of single- 
magnon excitation intensities in IPA-CuOU (symbols). Solid 
lines are simulations, as described in the text. 



main to decay into bunches of multiple magnons. For 
%^-/i c < h < ^ ti +/i c the magnons in IPA-CuCl 3 are 
thus stable. Their contribution to the dynamic struc- 
ture factor are sharp ^-function peaks. This restriction 
is removed for magnons inside the continuum: they are 
rendered unstable through decaying into pairs of lower- 
energy magnons. Their energy width (inverse lifetime) is 
dramatically increased. In IPA-CUCI3 this effect is par- 
ticularly strong, and for + h c < h < - h c there 
are no detectable single-magnon states. 

What drives the decay process is interactions between 
magnons. Interestingly, in an ideal symmetric spin lad- 
der the matrix elements between single-magnon and 2- 
magnon states are exactly zero, so the single mode sur- 
vives inside the 2-magnon continuum |20(. In particular, 
in an isotropic Haldane spin chain stable magnons are 
expected to persist all the way to the zone-boundary. 
We can speculate that the abrupt disappearance of the 
magnon branch in IPA-CUCI3 , which distinguishes these 
composite Haldane spin chains from actual ones, is due 
to a symmetry-breaking diagonal interaction J3. 

We finally want to point out the relevance of 
our results to the behavior of isostructural mixed 
(CH 3 )2CHNH 3 Cu(Cl a: Bri_ a: )3 materials In these 

compounds the spin gap decreases with increasing x, 
closes at x ~ 0.44 (the system becomes magnetically 
ordered), and then re-opens at x > 0.87. It was pro- 
posed that Br doping can change J\ from J\ < at 
x = 1 to Ji > at x = 12]. The phase diagram 
of such bond-disordered ladders has not been studied to 



date. However, on a naive mean- field level we can think 
of doping as a way to continuously tune J\. In this case, 
increasing x should lead the system through a gapless 
quantum-critical state, which for J3 = is exactly at 
Ji = 0. The reenter ant doping phase diagram is then 
easy to understand. Note that in the previously pro- 
posed alternating-chain model the gap remains open at 
Ji = 0. 

Contrary to what was previously thought, IPA- 
CUCI3 is a beautiful asymmetric spin ladder system, that 
for many purposes can be seen as a "composite" Haldane 
spin chain. Unlike an ideal Haldane spin chain though, it 
features a spectacular end-of-spectrum effect at a certain 
critical wave vector. Future neutron work will focus on a 
direct observation of excitation continua, effects of high 
magnetic field and Br-dilution. 
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